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Linear control theory is used to develop an improved lo- 
calized control scheme for spatially extended chaotic systems, 
which is applied to a Coupled Map Lattice as an example. The 
optimal arrangement of the control sites is shown to depend 
on the symmetry properties of the system, while their mini- 
mal density depends on the strength of noise in the system. 
The method is shown to work in any region of parameter space 
and requires a significantly smaller number of controllers com- 
pared to the method proposed earlier by Qu and Hu A 
nonlinear generalization of the method for a 1-d lattice is also 
presented. 

PACS numbers: 05.45. +b, 02.30.Wd 



Controlling chaos in high-dimensional systems S and 
spatiotemporal chaos especially is a very important prob- 
lem with numerous applications to turbulence || , insta- 
bilities in plasma [4], multi-mode lasers || and reaction- 
diffusion systems [6|. 

The present letter represents an effort to develop a 
general control algorithm for spatiotemporally chaotic 
systems using the methodology of linear control theory, 
which already proved to be fruitful Q] . Clarifying a num- 
ber of issues will have direct bearing on this. For in- 
stance, it is not clear how many parameters are required 
for successful control. If the control is applied locally, 
what is the minimal density of controllers and how should 
they be arranged to obtain optimal performance? What 
arc the limitations of the linear control scheme and how 
can they be overcome? 

Consider the Coupled Map Lattice (CML), originally 
introduced by Kaneko 0] : 
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with 



i — 1,2, ■■■,L and periodic boundary conditions, 
i.e. z\, L = z\ imposed. We also assume that the local 
map f(z, a) is a nonlinear function with parameter a, 
such that f(z*,a*) — z* . 
To be specific, we choose 



f(z) = az(l - z), 



(2) 



but emphasize, that all the major results hold indepen- 
dent of this choice. This CML has a homogeneous steady 
state z* — 1 — -, which is unstable for a > 3.0 and our 
goal is to stabilize it using a minimal number of con- 
trollers. 



The first attempt in this direction was undertaken by 
Hu and Qu jl| . The authors tried to stabilize the homo- 
geneous state by controlling an array of M periodically 
placed pinning sites ■ • ■ , %m} with appropriately cho- 
sen control u* 
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This however required a very dense array with distance 
between controllers L p = L/M < 3 in the physically 
interesting interval of parameters 3.57 < a < 4.0. 

The reason for this is the spatial periodicity of the 
pinnings. Since the system is spatially uniform, its eigen- 
modes are just Fourier modes and the pinning sites do 
not affect the modes whose nodes happen to lie at the 
pinnings, i.e. modes with periods equal to 2L p , 2L p /2, 
2L p /3, etc, provided those are integer. The control 
scheme worked only when all such modes were stable. 

It is however not necessary to destroy the periodicity 
completely to achieve control: that would complicate the 
analysis unnecessarily. Instead we add one more pinning 
site between each of the existing ones. Not all positions 
are good, but some do solve the problem — previously 
uncontrollable modes become controllable. 

In order to understand how the pinnings should be 
placed and see whether we achieve improved performance 
by introducing additional controllers, we have to use a 
few results of the linear control theory ||. We will start 
with linearizing eq. (|^) about the homogeneous steady 
state z* = (z*, ■ ■ ■ , z*) in both the state vector and con- 
trol to obtain the following standard equation 
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(4) 



where we denoted the displacement x = z — z*. If we 
define a = df(x*, a*)/dx, then the L x L Jacobian A is 
given by 
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and the L x M control matrix Bij = S(j — m)d(i — i m ) 
depends on how we place the pinning sites. 

If we use synchronous linear feedback u* = —Kx 1 , 
equation (Eh becomes 



x* +1 = (A-BK)^, 



(6) 
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and the solution x — can be made stable by a suitable 
choice of the feedback gain matrix K, if the controllabil- 
ity condition rank(C) — Lis satisfied. The controllability 
matrix C is defined via 

C = (B AB ■■■ A L ~ 1 B). (7) 

One can easily verify that the matrix B calculated 
for a periodic array of pinning sites does not satisfy the 
controllability condition and therefore the homogeneous 
steady state is not controllable. It can be stabilized if 
the weaker stabilizability condition is satisfied, i.e. all 
uncontrollable modes are stable. However this imposes 
excessive restrictions on the pinning density. 

The condition for stabilizability can be obtained from 
the spectrum of eigenvalues of the matrix (|^) 

7i = a(l-2e(l-cos(* i )) ) (8) 

where k\ = 0, fcj = fej+i = iri/L for i — 2, 4, 6, • • • and, 
for L-even, = 7r and a — 2 — a. Specifically, we need 
|(a-2)(l-2e(l-cos(7ri/i p ))| < 1 for all j = l,---,L-2, 
such that Lp/j is integer. Using this criterion one can 
obtain the relation between the minimum coupling, the 
distance between controllers and parameter a of the local 
chaotic map for a stabilizable system. For instance, j = 1 
yields 



2(a-2)(l-cos(7r/L p ))' v ' 

The results are presented in fig. [l| It can be easily ver- 
ified that they coincide with the numerically obtained 
results of Hu and Qu. 

It is possible however to extend the limits of the control 
scheme quite substantially by making the system control- 
lable as opposed to stabilizable. This is easily achieved 
by choosing a different matrix B, i.e. placing the pin- 
ning sites differently. Doing so will enable us to control 
the system anywhere in the parameter space at the same 
time using a smaller density of controllers. 

First one has to determine the dimensionality of the 
matrix B, in other words determine the minimal num- 
ber of parameters required to control the CML (|l|) of an 
arbitrary length. It can be shown fic|| that the minimal 
number of parameters required to control a system with 
degenerate Jacobian is equal to the greatest multiplicity 
of its eigenvalues. 

Since the system under consideration has parity sym- 
metry, the eigenvalues (||) of its Jacobian are in fact dou- 
bly degenerate, so the minimal number of control param- 
eters yielding a controllable system in our case is two, 
meaning at least two pinning sites are required. One can 
easily verify that the controllability condition for anLx2 
matrix 

Bn = 5(j - l)5(i - ii) + 5(j - 2)5{i - i 2 ) (10) 

is indeed satisfied for a number of arrangements {ii,^}- 
The restrictions on the mutual arrangement of the con- 
trollers are again given by the condition of controllability: 
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FIG. 1. Periodic array of single pinning sites: critical 
coupling e cr as a function of parameter a. The dots represent 
the numerical results from figure 2 of [1], with e rescaled by 
a factor of two to make it compatible with our definition. 

L should not be a multiple of |i2 — i\ | , otherwise the mode 
with the period 2\%2 — i\ \ becomes uncontrollable. 

The next step in the algorithm is to determine the 
feedback gain K. Pole placement techniques based on 
Ackermann's method [fLlJ are inapplicable to the problem 
of controlling spatially extended systems because they 
are numerically unstable |]l2|| and break down rapidly for 
problems of order greater than 10. 

Instead we use the method of the linear-quadratic (LQ) 
control theory Q , applicable to the unstable periodic tra- 
jectories as well as fixed points. This method is not only 
numerically stable, but also allows one to optimize the 
control algorithm to increase convergence speed, and at 
the same time minimize the strength of control. As we 
will see below, decreasing control enlarges the basin of at- 
traction, which has very important consequences for the 
time to achieve control (capture the chaotic trajectory). 
The optimal solution is obtained by minimizing the cost 
functional 

oo 

V(x°) = ^(x it Qx* + u*W), (11) 

n=0 

where Q and R are the weight matrices that can be cho- 
sen as any positive-definite square matrices. 
The minimum of ( fj"l| ) is reached when 

K = (R + B^PBY l B^PA 1 (12) 

where P is the solution to the discrete-time algebraic 
Ricatti equation 

P = (Q + A" 1 PA) - A^P^B(R + B^PB^B^PA. (13) 

Numerical simulations show that the CML (Q,§) can 
indeed by stabilized by this linear control scheme in a 
wide range of parameters a and e. The solution for K is 
presented in figure | for a = 4.0, e = 0.33 and L — 8 with 
Q = hxs and R = 12x2- The steady homogeneous state 
z* = 0.75 has 3 unstable and 5 stable directions and we 
use 2 pinning sites to control it. 
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FIG. 2. Feedback matrix: feedback gain vectors K\j and 
K23 for left and right controller placed at the sides of the 
lattice (i± = 1,«2 = 8) as functions of the lattice site for 
a = 4.0 and e = 0.33. 



The contribution —K m iX™ from the sites i far away 
from the pinning site i m is larger, as one would expect: 
since the feedback is applied indirectly through coupling 
to the neighbors, the perturbation introduced by the con- 
trollers decays with increasing distance to the pinning 
sites. 

Noise limits our ability to locally control arbitrarily 
large systems with local interactions. We will use a sim- 
ple illustrative approach to see the effect of noise on the 
control scheme. The rank of the matrix is given by 
the number of its nonzero singular values. The singu- 
lar values of the controllability matrix (Q) scale roughly 
as si ~ |7i| ( , where 71 is the largest eigenvalue of the 
Jacobian A 
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e < 0.5, 
e > 0.5. 



(14) 



Assuming that there is an uncertainty in the calculations 
(due to the uncertainty in the state vector, parameter 
vector or just numerical roundoff errors) of relative mag- 
nitude er, we can say that the rank of the controllabil- 
ity matrix can be reliably determined to be equal to the 
length of the lattice if sq/sl > o. This gives us the the- 
oretical bounds on the size of the controllable system in 
the presence of noise: 



L« = 



logO) 



(15) 



On the other hand, the perturbation 8xi introduced by 
the controller i affects the dynamics of the remote site j 
after propagating a distance A = \i — j\ in time r = A, 
decaying by a factor of e per iteration, while the noise at 
site j increases roughly by a factor of 71 per iteration. 
We therefore need Sxie A > <t|7i| t . Since the maximum 
distance A to the closest controller is L/2 and Sxi ~ 1, 
we get another bound, complementing dl5l) 



21og(<r) 



log(e) - A,, 



(16) 
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FIG. 3. The largest length of the lattice which can be 
controlled with two pinning sites: theoretical estimates (solid 
lines) and numerical results (dots) obtained with the uni- 
formly distributed noise of amplitude a — 10~ 14 as functions 
of coupling e for a = 4.0. 



Similar constraints were obtained by Aranson et. al for 
the lattices with asymmetric coupling (cf. equation (15) 
of the ref. @). 

The maximal length of the system, that can be stabi- 
lized by the LQ method with two pinning sites placed 
next to each other is obtained numerically by choos- 
ing the initial condition very close to the fixed point 



(|x°| 



l-L/2 



) and letting the system evolve under 



control (TL2f) calculated for Q = Ilxl and R = 12x2- 
This length is quite large even in the presence of noise 
(fig. ||) and agrees with the theoretical bounds (pl| |l^) 
rather well for such a crude estimate. 

The problem of controlling a large 1-dimensional sys- 
tem with the length L > L max (cr) exceeding the max- 
imum allowed for a given noise level can be easily re- 
duced to the problem of controlling a number of smaller 
systems with the length L p < L max (a). We partition 
the entire lattice {z{, 
mains {z\ m _ 1)Lp+v ■■ ■ 
ray of pinning sites 



1ml 



z l L } into M = L/Lp subdo- 
lLp }, and control it with an ar- 
= (m - l)L p + 1, i m2 = mLp, 



m = 1, • • • M positioned periodically at the boundaries 
of these subdomains. 

The stabilization can be achieved by choosing 

Lp Lp 

+ Jj0(<fot - |ff( m -i)L„+»l) ^ Kli x \m-l)L p +i 
i=l i=l 

U im2~ F( Z l m2 ~l-< Z i m2 l Z i ml ) — F\ Z i m2 -ll Z l ra 2' Z im2+l) 

Lp Lp 

+ ]J 9 ( 6x i - \ X \m-l)L p +i\) J2 K ^ X \m-l)L p +i, ( 17 ) 
i=l i=l 

where 6{x) is a step-function. 

This arran gem ent effectively carries two functions. We 
use control (0) to (nonlinearly) decouple the subdo- 
mains, simultaneously imposing periodic boundary con- 
dition for each subdomain (the first two terms) to make 
the system controllable. Then we stabilize each subdo- 
main asynchronously by applying a linear (in deviation 
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FIG. 4. Stabilizing uniform steady state: a large lattice 
(L — 128) is controlled by an array of double pinning sites, 
placed at the boundaries of subdomains with length L p — 8. 
The state of the system was plotted at each 10000-th step. 



x\ = z\ — z*) feedback (the last term), inside the neigh- 
borhood of the fixed point determined by bxi . The linear 
approximation (^) is only valid if 



To summarize, we have shown that the restrictions on 
the minimal density of periodically placed single pinning 
sites obtained by Qu and Hu H as a result of numerical 
simulations can in fact be obtained analytically from the 
stabilizability condition. 

The efficiency of the control scheme can be improved 
significantly if one uses double pinnings instead of single 
ones. The homogeneous steady state becomes control- 
lable for any values of the control parameters and the 
minimal density of pinning sites is reduced substantially. 
It is shown that the maximal distance between the pin- 
nings depends on the strength of noise in the system and 
can be estimated analytically. 

The appropriately chosen (using the LQ technique) 
feedback can decrease the capture time for the chaotic 
trajectory by enlarging the capture region. The introduc- 
tion of nonlinearity into the control scheme can decrease 
this time even more significantly by effectively decoupling 
the large lattice into a number of smaller subdomains. 

The authors thank Prof. J.C. Doyle for many fruitful 
discussions. This work was partially supported by the 
NSF through grant no. DMR-9013984. H.G.S. thanks C. 
Koch for the kind hospitality extended to him at Caltech 
and the Volkswagen Foundation for financial support. 
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and therefore strong feedback significantly decreases the 
size of the capture region, which makes the capture 
time vary large. Minimizing the capture time can be 
achieved by minimizing the feedback strength using the 
LQ method (|l|||). 

We demonstrate this approach by stabilizing the homo- 
geneous stationary state of the CML defined by equation 
(|,|) with a = 4.0, e = 0.33. L = 128 sites were di- 
vided into M = 16 subdomains of length L p = 8, each 
controlled by two pinning sites. The results presented in 
fig. | show the evolution of the system from the initial 
condition chosen to be a collection of random numbers 
in the interval [0, 1]. 

Eqs. ( |i~5| , [l6| ) now give the minimal density of pinning 
sites that yields the controllable fixed point solution. It 
is indeed seen to be much lower than that given by (|^), 
e.g. 2/L p = 1/20 (1/17 from the numerics, see fig. ||) as 
opposed to 1/Lp = 1/2 for the choice a = 4.0, e = 0.4 
and the precision of calculations given by a = 10~ 14 . 

Although the resulting control scheme becomes non- 
linear (and therefore requires full knowledge of the evo- 
lution equations), it has the additional benefit, that the 
capture time is determined by the length L p <C L and 
is typically many orders of magnitude smaller than that 
obtained for the linear control scheme (obtained by lin- 
earizing (|l7j)), which only requires the Jacobian to be 
known. In fact our computational resources were insuffi- 
cient to observe even a single capture for L > 40 with the 
linearized control. Generalizing this nonlinear approach 
to higher-dimensional systems remains a challenge. 
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